21-cm Background Anisotropies Can Discern Primordial Non-Gaussianity 
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The non-Gaussianity of initial perturbations provides information on the mechanism that gen- 
erated primordial density fluctuations. We find that 21-cm background anisotropies due to inho- 
mogeneous neutral hydrogen distribution prior to reionization captures information on primordial 
non-Gaussianity better than a high-resolution cosmic microwave background anisotropy map. An 
all-sky 21-cm experiment over the frequency range from 14 MHz to 40 MHz with angular information 
out to a multipole of 10 5 can limit the primordial non-Gaussianity parameter /nl <j 0.01. 
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Introduction — The cosmic 21-cm background involving 
spin-flip line emission or absorption of neutral hydro- 
gen contains unique signatures on how the neutral gas 
evolved from last scattering at z ~ 1100 to complete 
reionization at z ~ 10 [l[. Subsequent to recombination, 
the temperature of neutral gas is coupled to that of the 
CMB. At redshifts below ~ 200 the gas cools adiabati- 
cally, its temperature drops below that of the CMB, and 
neutral hydrogen resonantly absorbs CMB flux through 
the spin-flip transition 0, 0] • The inhomogeneous neu- 
tral hydrogen density distribution generates anisotropies 
in the brightness temperature measured relative to the 
blackbody CMB 0. The large cosmological information 
content in 21-cm background is well understood in the 
literature 0, [||. If the primordial fluctuations are 
non-Gaussian, then the 21-cm anisotropies will naturally 
contain a signature associated with that non-Gaussianity. 

In this Letter, we show that the angular bispec- 
trum of 21-cm anisotropies contains a measurable non- 
Gaussianity from primordial fluctuations, if the scale- 
independent quadratic corrections to the gravitational 
potential captured by /nl has a value above ~ 0.01. This 
is substantial given that even a perfect CMB experiment 
limited by cosmic variance alone can only restrict /nl > 3 
while there is no significant improvement when using 
low redshift large-scale structure [8| . In comparison, the 
expected non-Gaussianity under standard inflation is of 
order \n s — 1|, and with the scalar spectral index n s ~ 0.98 
[lH . this non-Gaussianity is well below unity j9|]. The 
primordial non-Gaussianity parameter /nl, however, has 
a correction associated with evolution of second and 
higher-order perturbations after inflation [l~oj | . For stan- 
dard inflation / NL = -5/12(n s - 1) + 5/6 + 3/10/(fc) 
where the last term is momentum dependent. In this 
case, a primordial non-Gaussianity is always present and 
/nl is at the level of a few tenths. 

Thus, compared to all other probes of inflationary pa- 
rameters, such as the spatial curvature and amplitudes 
and spectral indices of scalar and tensor spectra, /nl is 
the only parameter for which we do not yet have an ob- 
servational probe to reach the simple expectation. Here, 
we establish 21-cm background as a potentially useful 



probe to study /nl values below one. When compared 
to CMB, 21-cm background has two distinct advantages: 
(1) the ability to probe multiple redshifts based on fre- 
quency selection and (2) the lack of a damping tail in the 
21-cm anisotropy spectrum, similar to early damping in 
CMB spectrum at a multipole around 2000. 

Here, we also consider the confusion from other non- 
Gaussian signals in the 21-cm sky. We introduce the 
two-to-one correlator statistic [l2| that can be optimized 
to detect the primordial bispectrum by appropriately fil- 
tering 21-cm anisotropy data. This statistic has already 
been used to limit /nl with CMB data [l3( . The statistic 
is optimized with a filter to extract most information on 
the amplitude of the primordial non-Gaussianity, but re- 
quires prior knowledge on the configuration dependence 
of the primordial 21-cm bispectrum. The amplitude of 
non-Gaussianity remains a free parameter to be deter- 
mined from the data. The Letter is organized as follow- 
ing: w e first discuss the bispectrum in 21-cm anisotropics 
associated with primordial perturbations resulting from 
quadratic corrections to the primordial potential. We 
discuss ways to measure this bispectrum in the presence 
of other non-Gaussian signals and determine the extent 
to which /nl can be measured from 21-cm background 
data. 

21-cm Signal — The 21-cm anisotropies are observed as 
a change in the intensity of the CMB due to line emission 
or absorption at an observed frequency v: 
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where T$ is the spin temperature of the neutral gas, 
z is the redshift corresponding to the frequency of ob- 
servation (1 + z = v 2 \/v, with v 2 \ = 1420 MHz) and 
Tqmb = 2.73(1+2:)^ is the CMB temperature at redshift 
z. The optical depth, r, in the hyperfine transition [2|, 
when accounted for density and velocity perturbations of 
a patch in the neutral hydrogen distribution 0,13, li | , is 
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where Aiq is the spontaneous emission coefficient for the 
transition (2.85 x 10~ 15 s^ 1 ), nni is the neutral hydrogen 
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FIG. 1: (a) The 21-cm anisotropy power spectrum 
(y/l 2 Ci /2ir; solid line) and non-Gaussian squared 21-cm 
anisotropy — 21-cm anisotropy power spectrum ([l 3 Xi/2n] 1 ^ 3 ; 
see text for details) with the optimal filter applied (see, equa- 
tion 12). The dashed line is the signal related to the pri- 
mordial bispectrum with /nl = 1, while the dotted line is 
the noise bias Xf rSlV from perturbation theory bispectrum un- 
der non-linear gravitational clustering. The dot-dashed line 
shows the total noise spectrum, [1 3 Ni/2ty] 1 ^ 3 , associated with 
the measurement of Xf rlm . These spectra are calculated at 
z — 100. (b) The cumulative signal-to-noise ratio for a de- 
tection of the non-Gaussian spectrum Xi as a function of the 
multipole I at z = 100 (solid line) and at z — 30 (dashed line) . 
These estimates assume no detector noise with measurements 
limited by cosmic variance alone and with / a k y = 1. 



density, Sh — (hh — ^h)/^h is the inhomogeneity in 
the density, v is the peculiar velocity of the neutral gas 
distribution, r is the comoving radial distance, and H(z) 
is the expansion rate at a redshift of z. For simplicity, 
we have dropped the dependences in location n. The 
fluctuations in the CMB brightness temperature is 
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where S(z) describes the coupling between fluctuations in 
the spin temperature and the neutral density distribution 
Q. In above % ps 26.7mKvTTi. 

In terms of non-Gaussian measurements, we discuss 
the angular averaged bispectrum of 21-cm anisotropics 
constructed with spherical harmonic moments of the 21- 
cm background defined as 
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For simplicity in notation, we will drop the explicit de- 
pendence on the frequency, but it should be assumed that 



21-cm background measurements can be constructed as a 
function of frequency, and thus as a function of redshift, 
with the width in frequency space primarily limited by 
the bandwidth of a radio interferometer. 

To generate a non-Gaussianity in the 21-cm back- 
ground, we assume quadratic corrections to the primor- 
dial potential fluctuations such that $(x) = $l(x) + 
/nl($|(x)-(<&|(x))) wnen $l( x ) is the linear and Gaus- 
sian perturbation and /nl is the coupling constant, which 
may or may not be scale dependent. Under this descrip- 
tion, existing WMAP data limits -30 < /nl < 74 at the 
la confidence while with an ideal CMB experiment 
fundamentally limited by cosmic variance one can con- 
strain |/nl| < 3 [7J. As discussed, this is well above the 
standard expectations for /nl even after accounting for 
second-order evolution during horizon exit and re-entry 

Using the resulting three-dimensional bispectrum for 
potential fluctuations, and noting that 21-cm observa- 
tions trace the density field, we can write the expected 
angular bispectrum of 21-cm background as 
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with Ji{x) = ji(x) — j'i(x). In equation ([BJ, D{r) is 
the growth function of density perturbations, M{k) — 
—3k 2 T(k)/5Q m HQ maps the primordial potential fluc- 
tuations to that of the density field, and W v {r) is the 
window function of observations related to the central 
frequency and the bandwidth. In equation permuta- 
tions account for a total of three terms with the replace- 
ment of I2 — > ^3 and l\. For reference, the angular power 
spectrum of 21-cm anisotropies follows from above as 
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While we have explicitly written out the bispectrum re- 
lated to primordial non-Gaussianity, 21-cm anisotropics 
will contain additional signatures of non-Gaussianity. 
Among these is the 21-cm bispectrum resulting from the 
non-linear gravitational evolution of density perturba- 
tions with a three-dimensional bispectrum of the form 
Bs(k!,k 2 ,k 3 ) = F(ki,k^)Pss(ki)Pss(k2 ) + Perm., where 
F(ki,k 2 ) is the mode-coupling kernel [15(. We denote 
the resulting angular bispectrum of 21-cm anisotropics 
sourced by the non-linear density field as Bf a ; when 
the above three dimensional bispectrum is substituted to 
the formalism outline with equation (5) Moreover, 
once the universe begins to reionize, 21-cm anisotropics 
are modulated by fluctuations in the neutral fraction, Xh, 
in addition to density perturbations with r cx n^x^l + 
<5h)(1 + <5 X )- Though small in amplitude, the second- 
order corrections to the brightness temperature lead to 
a non-Gaussian signal 17). As we are concentrating on 
redshifts prior to reionization, this secondary bispectrum 
can be safely ignored. 

For a measurement of /nl we consider a simple statis- 
tic that captures non-Gaussian information once an ex- 
pected configuration for the 21-cm bispectrum is spec- 
ified. The statistic involves the power s pec trum con- 
structed with the filtered 2-to-l correlator [l2l.[l3j]: 
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where af m = J dnT fc 2 (n)Y, m *(n), where f 6 (n) is the 21- 
cm brightness temperature with the filter described be- 
low applied in multipole space. The measurable statistic 
can be written as 
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where (i) represents the 21-cm bispectrum generated by 
either primordial non-Gaussianities or non-linear gravi- 
tational evolution. Since the bispectrum is defined by a 
triangle in multipole space with lengths of sides (Z, Zi, Z2), 
the Xi statistic essentially captures information from all 
triangular configurations of the bispectrum with one of 
the sides fixed at length I. The filter function w(l\,l2\l) 
in multipole space, however, is applied so that one can 
maximize the information related to one form of non- 
Gaussianity with a specified configuration dependence. 
Here, this allows a mechanism to separate the primordial 
non-Gaussianity from other non-Gaussian signals since 
we have prior theoretical expectation on the the expected 
configuration dependence [lOj]. 

The signal-to-noise ratio for detecting Xf rlvn with an 
experiment covering / s k y fraction of the sky is S/N = 



/sky E*(^? nm )7W ot ) 2 > where ^i tot is the total noise 
contribution which involves both the noise bias from 
Af lav and the Gaussian variance in the proposed statis- 
tic: 
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Here, C\ ot represents all contributions to the 21-cm 
anisotropy power spectrum such that C\ ot = C ; 21 ~ cm + 
^rnoise w \ ien Qioise j g instrumental noise contribu- 
tion. Note that we are treating the confusion, X ; grav , in 
the same manner foreground signals are treated in CMB 
studies. 

As suggested above, to maximize the signal-to-noise 



ratio for a detection of A ; prim , following Ref. |12| we found 
the filter to be of the form 
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where &f;™ 2 = A^i 1 j 2 . The filter is normalized 

such that Y~\ 1 wP rl , m ,, = 1. The filter is written in terms 

^^'1*2 tl,t2|t 

of Bf\" n l2 , but due to normalization of the filter, is inde- 
pendent of /nl, thus no prior knowledge on the ampli- 
tude of the non-Gaussianity is required. What is neces- 
sary, however, is the expected configuration dependence, 
which can be specified a priori under an assumed theo- 
retical model for primordial non-Gaussianity. 

Since the filter prescribes the expected bispectrum con- 
figuration, it also acts as a way to reduce the ampli- 
tude of confusions. For example, with the filter ap- 
plied, Xf la,v oc Yli i 2 Bf ; a \ Bfi™ 2 an d as modes do not 
align exactly, there is a significant cancellation. We find 
that the resulting noise spectrum for a measurement of 
Xf mm is dominated by the Gaussian variance captured 
by the second term of equation (fTT|) and not the con- 
fusion resulting from Xf Tav . This statement is indepen- 
dent of /nl ■ If stated differently, when properly filtered 
to search for the primordial non-Gaussianity, the main 
confusion for detecting primordial signal is not the non- 
Gaussianity generated by non-linear perturbations but 
rather the Gaussian covariance associated with the sta- 
tistical measurement of A ; prim . 

Finally, while we have not explored here as the Gaus- 
sian variance dominates, it may be possible to further 
optimize the filter such that the additional non-Gaussian 
confusions cancel out exactly. In practice, the measure- 
ment of X prim is likely to be further confused by the 
non-Gaussianity of foregrounds. Unfortunately, little is 
known about the expected level of the foreground inten- 
sity in the frequency range of interest. Techniques have 
been suggested and discussed to remove foregrounds be- 
low the detector noise levels flil . [l9| and the filtering 
process we have outlined will further reduce the confu- 
sion from the remaining residual foregrounds. This is 
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clearly a topic for further study once data become avail- 
able with first-generation interferometers 1 . 
Results — Fig. 1(a) summarizes the power-spectrum of 
21-cm anisotropies generated by the neutral hydrogen 
distribution at a redshift of 100. Here, we also plot Xf llm , 
X, grav , and Ni for the same redshift with the optimal filter 
applied. As shown, Ni > Xf rav , suggesting that the noise 
term is dominated by the Gaussian variance (second term 
in equation 11). Note that we have estimated JVj under 
the assumption that observations are limited only by the 
cosmic variance. This allows us to establish the poten- 
tially achievable limit and compare directly with cosmic 
variance limit of /nl with CMB. 

In Fig. 1(b), we summarize the estimate related to 
signal-to-noise ratio for a detection of X ; prim as a func- 
tion of /. The typical signal-to-noise ratio, when mea- 
surements are out to a multipole of 10 5 , is at the level 
of ~ 13/nl, with observations centered at a redshift of 
100, to ~ 5/nl, at a redshift of 30. If one restricts the 
calculation to certain configurations, such as equilateral 
triangles of the bispectrum, then the cumulative signal- 
to-noise ratio is ~ 0.3/nl [20(] . Out to a multipole of 10 5 , 
however, there is a total of ~ 10 15 modes in the whole 
bispectrum, whereas with a specific configuration, one 
is only restricting to at most 10 5 modes. Our statistic 
X[ captures all this information and is optimally filtered 
so as to maximize information on the primordial non- 
Gaussianity. When using all the information, though 
there are 10 15 modes, the filtering selection removes a 
large number of modes when separating the primordial 
signal from the secondary non-Gaussian signal. 

Furthermore, 21-cm observations lead to measure- 
ments at multiple redshifts. One cannot, however, make 
arbitrarily small bandwidths to improve the detection 
since at scales below a few Mpc, anisotropies in one red- 
shift bin will be correlated with those in adjacent bins 
fl9j |. With a bandwidth of 1 MHz, corresponding to a 
radial distance of 30 Mpc at a redshift of ~ 50, over 
the frequency interval between 14 MHz (z ~ 100) to 45 
MHz (z ~ 30), we find a cumulative signal-to- noise ratio, 
S/N = VE/[S/N(z)] 2 , of ~ 10 2 / NLv <fe out to l max of 
10 5 . In return, one can potentially probe /nl values as 
low as 10 -2 . With first-generation instruments, we would 
at most survey 1% of the sky. Assuming instrument noise 
is dominating at multipoles above 10 4 between 30 MHz 
and 60 MHz at 1 MHz bandwidths, we find a signal-to- 
noise ratio of 0.5/nl\/ /sky /0. 01, which is equivalent to 
the ultimate limit from CMB. 

The lowest /nl value one can reach with 21 cm data is 
two orders of magnitude better than CMB [3] . This is the 
main result of the paper with important implications for 
future cosmological studies. While it may be challenging 

1 http://www.lofar.org; http://www.haystack.mit.edu/arrays/MWA 



for first-generation interferometers, low frequency 21-cm 
observations should eventually establish the primordial 
non-Gaussianity. The limitation to achieve this detec- 
tion is not a fundamental one, such as cosmic variance, 
but rather that involving adequate technology for high- 
sensitive low-frequency radio observations. In the same 
manner technology has been steadily improving over the 
last two decades to allow high precision CMB measure- 
ments, it is very likely that experiments will progress 
rapidly allowing precise 21-cm studies. While we lacked 
a unique probe to reach /nl values below 0.1, our cal- 
culations suggest that this is no longer the case once we 
open the low-frequency radio sky for observations. 
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